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ABSTRACT 

We report on the discovery of a high- velocity narrow absorption line outflow in 
the redshift 2.3 quasar J212329. 46-005052. 9. Five distinct outflow systems are detected 
with velocity shifts from -9710 to -14,050 km s _1 and C iv AA1548,1551 line widths 
of FWHM w 62 to 164 km s" 1 . This outflow is remarkable for having high speeds and 
a degree of ionization similar to broad absorption line (BAL) flows, but line widths 
roughly 100 times narrower than BALs and no apparent X-ray absorption. This is also, 
to our knowledge, the highest-velocity narrow absorption line system confirmed to be 
in a quasar outflow by all three indicators of line variability, smooth super-thermal 
line profiles and doublet ratios that require partial covering of the quasar continuum 
source. All five systems have stronger absorption in O vi AA1032,1038 than C iv with no 
lower ionization metal lines detected. Their line variabilities also appear coordinated, 
with each system showing larger changes in C iv than O vi and line strength variations 
accompanied by nearly commensurate changes in the absorber covering fractions. The 
metallicity is approximately twice solar. 

These data require five distinct outflow structures with similar kinematics, physi- 
cal conditions and characteristic sizes of order 0.01-0.02 pc (based on partial covering). 
The coordinated line variations, occurring on time scales ^0.63 yr (quasar frame), are 
best explained by global changes in the outflow ionization caused by changes in the 
quasar's ionizing flux. An upper limit on the acceleration, <3 km s _1 yr -1 , is con- 
sistent with blobs of gas that are gravitationally unbound and coasting freely >5 pc 
from the central black hole. Additional constraints from the variability time indicate 
that the full range of plausible distances is 5 < R < 1100 pc. However, if these small 
absorbing structures were created in the inner flow, they should be near the ~5 pc 
minimum radius because they can travel just a few pc before dissipating (without ex- 
ternal confinement). An apparent double line- lock in ClV suggests that the flow was 
radiatively accelerated and its present trajectory is within ^16° of the radial (line-of- 
sight) direction. The absence of strong X-ray absorption shows that radiative shielding 
in the far-UV and X-rays is not needed to maintain moderate BAL-like ionizations 
and therefore, apparently, it is not needed to facilitate the radiative acceleration to 
high speeds. We argue that the ionization is moderated, instead, by high gas densities 
in small outflow sub-structures. Finally, we estimate that the kinetic energy yield from 
this outflow is at least two orders of magnitude too low to be important for feedback 
to the host galaxy's evolution. 

Key words: galaxies: active — quasars: general — quasars: absorption lines - 
quasars: individual: J212329.46-005052.9 



2 F. Hamann et al. 



1 INTRODUCTION 

High-velocity outflows from quasars appear to be a nat- 
ural part of the accretion process. They are detected 
most conspicuously via broad absorption lines (BALs) in 
quasar sp ectra, with velocity w idths >2000 km s _1 (by 
definition, Wc vmann et al.l Il99ll ) and flow speeds from a 
few thousand to roughly 10,000-30,000 km s _1 (see also 
iKorista et all 1 19931 ; iTrump et all l2006h . However, quasar 
outflows can also produce narrow absorption lines (NALs), 
with velocity widths less than a few hundred km s , 
as well as mini-BALs, which are loosely defined as inter- 
medi a te between NALs a n d BALs (e.g., Hamann fc Sabral 
12004 iNestor et all l20Qg| ; IRodriguez Hidalgo et all |2010| ). 
These narrower outflow lines are more common than 
BALs (IRodriguez Hidalgo et al.1 l20ld : IGibson et al.l 
Nestor et all 120081 ; IWild et all 120081 ; iMisawa et all ' 
Richards! 120011 ). Overall it is estimated that >50% of 
optically-selected quasars have at least one type of 
outfl ow absorption line in the i r rest-frame UV spec- 
tra ( Ganguly fc Brothertonl 120081 ; IRodriguez Hidalgo et al.l 
I201C ) . The flows themselves are probably present in all 



2009; 



2007; 



quasars if, as expected, the outflowing gas covers just a frac- 
tion of the sky as see n from the central emission source (e.g., 
lHamann et~atlll993h . 

Quasar outflows have gained attention recently as 
a mechanism that can physically couple quasars to the 
evolution of their host galaxies. In particular, some 
models of galaxy evolution invoke the kinetic energy 
"feedback" from an accreting super-massive black hole 
(SMBH) to regulate both the star formation in the 



tral SMBH dKauffmann & Haehnerd 2000; 


Granato et al.l 


l2004l;lDi Matteo et al.ll2005l; Hopkins & Elvisl 


2010, and refs. 



therein). This type of coupling via feedback could pro- 
vide a natural explanation for the observed mass corre- 
lation between SMBHs and their host galaxy spheroids 
ilTremaine et al.ll2002l ; iMarconi fc Huntll2003l ; lHaring fc Rixl 



1 2004 ) . The kinetic energy luminosity, Lk , needed for signif- 
icant feedback from a quasar outflow is believed to be just 
a few percent of t he quasar bolometric (photon) luminosity, 
e.g., Lk/L ~ 5% dScannapieco fc Obll2004l ; lDi Matteo et all 



120051 : IProchaska fc Hennawil boOO ) . 



Unfortunately, many aspects of the energetics and phys- 
ical properties of quasar outflows remain poorly understood. 
Computational models aimed mainly at BALs attribute the 
flows to gas lifted off of the accretion disk and driven out- 
ward to high spee ds by radiati o n pressure or magneto- 



centr if ugal forces dMurrav et all 19951; iMurrav fc Chiang! 



1 19971 ; iProga fc Kallmanl |2004| ; lEverettl l2005h . Radiative 
forces are expected to dominate in quasars, where the bolo- 
metric luminosi ties are a subs tantial fraction of the Edding- 
ton luminosity l|Everettll2005T ). 

One essential requirement for radiative driving is that 
the outflow is not too highly ionized, so it can maintain 
significant line and bound-free continuum opacities across 
the main acceleration region. This is not trivial for flows 
launched from the inner accretion di sk, near the quasar's 
intense source of ionizing radiation. IMurrav et ail (|l995l ) 
proposed that the outflow ionization is moderated by a 
high column density of roughly stationary material that 
serves as a radiative shield at the base of the flow (also 



IMurrav fc Chiangll 19971 ). The shielding medium is itself too 
ionized and too transparent to be driven radiatively to 
high speeds, but it might have enough continuous opac- 
ity in the far-UV and X-rays to lower the ionization and 
thereby facilitate the acceleration of the gas directly behind 
it (i.e., farther from the emission source). This scenario is 
supported by observations showing that BAL quasars are 
heavily obscur ed in X-rays, consisten t with strong radia- 
tive shielding (iGreen fc Mathuij Il996l ; iMathur et all |2000| ; 
iGallagher et al.ll2002l . 2006). However, other observations in- 
dicate that quasars with NAL or min i-BAL outflow s have 
substantially less X-ray absorption l|Chartas et ail 120091 ; 
iGibson et all 120091 ; IMisawa et ail 120081 ). This disparity in 
the X-ray results, between BALs on the one hand and NALs 
and mini-BALs on the other, presents an important chal- 
lenge to the paradigm of radiative acceleration behind an 
X-ray/far-UV shield (see also §4.5 below). 

A central problem is our poor understanding of the 
physical relationships between BAL, NAL and mini-BAL 
outflows. One possibility is that the different line types 
are simply different manifestations of a single outflow phe- 
nomenon viewed at different angles. For example, it is of- 
ten supposed that BALs form in the main body of the out- 
flow near the accretion disk plane, while NALs and mini- 
BALs form along sightlines that skim the edges of the BAL 
flow at higher latitu des above the disk l|Gangulv et al.ll200ll ; 
IChartas et al.l [2009) . There could also be evolutionary ef- 
fects. For example, NALs and mini-BALs might represent 
the tentative beginni ng or end stages of a more powerful 
BAL outflow phase (|Hamann et alj|2008l ). There is some 
evidence for outflow evolution in that a particular variety 
of low-ionization BALs (the so-called FeLoBALs) is found 
preferentially in dusty young host galaxies with high star 
formation rates l|Farrah et all 2007). Variability studies have 
shown further that BALs and mini-BALs in quasar spectra 
can appear, disappear or swap identities (where a BAL be- 
comes a mini-BAL or vice versa) on time scales of months 
to years in t he quasar rest frame llHamann et al. 20081; 



Gibson et all 120081; iLeighlv et al.1 12009]; ICapellupo et all 
20101; IRodriguez Hidalgo et al.1 |2010| ; IGibson et all l20l6T T 
The true relationship between outflow NALs, BALs and 
mini-BALs probably involves an amalgam of orientation and 
time-dependent effects. 

To build a more complete picture of quasar outflows, 
we need better observational constraints on the physi- 
cal properties of each outflow type. NAL outflows present 
a unique challenge because they are difficult to identify. 
Most narrow absorption lines in quasar spectra form in 
cosmologically intervening gas that has no relationship 
to the background quasar. Statistical studies that exam- 
ine the velocity d i stributions of the narrow metal lines 
l|Nestor et al.1 120081 ; I Wild et al.l l200Sf ). or correlations be- 
tween the incidence of th ese lines and t he properties of 
the background quasars l|Richardsl l200ll ), indicate that 
significant fractions of the narrow metal-line systems do 
arise in quasar outflows. However, to identify individ- 
ual NAL outflow systems we need more information. The 
most commonly used indicators of an outflow origin are 
line variability, resolved absorption profiles that are sig- 
nificantly broad and smooth compared to thermal line 
widths, or line strength ratios in multiplets that reveal par- 
tial line-of-sight covering of the background light source 
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IIHamann et al. 1997 ; Barlow fc SargentJI 19971; iMisawa et al.l 



2007; lHamann fc Simo n 2010; lSimon fc Hamanrj|2010f) . 

Here we discuss a complex of five outflow NAL systems 
at velocities from —9710 to —14,050 km s" 1 in the redshift 
~2.3 quasar J212329.46-005052.9 (hereafter J2123-0050). To 
our knowledge, these are the highest velocity NALs whose 
outflow nature is confirmed by all three indicators of vari- 
ability, partial covering and smooth super-thermal line pro- 
files. We selected this quasar from the Sloan Digital Sky Sur- 
vey (SDSS) as part of a larger study of intervening, metal- 
strong d amped Lycv (DLA) and sub-DLA absorption line 
systems jH crbcrt-F ort et all 2006f). Follow-up observations 
jKaplan fc et al.ll2010l ; iMalec et al.ll2010l ; iMilutinovic et all 
|2010| ) revealed the variability and outflow origin of the NAL 
systems discussed here. In the sections below we describe 
the observations (§2), present measurements and analysis of 
the outflow lines (§3), and discuss the interpretation and 
broader implications of our results (§4). We provide a de- 
tailed general summary in §5. Throughout this paper, we 



Table 1. Observation Summary 



adopt a cosmology with H 
and Qa = 0.7. 



70 km s _i Mpc, ftj, 



0.3 



2 OBSERVATIONS & OTHER DATA 

Table 1 summarizes the observations used in this study. 
Columns 1 and 2 give the fractional years and calendar 
dates of the observations. Columns 3-5 list the telescope, 
spectral resolution and wavelength coverage, AA, respec- 
tively. We obtained the 2002.68 spectrum from the SDSS 
archives fully reduced and calibrated. The fluxes provided 
by the S DSS are expected to be accurate to within a few 
percent (|Adelman-McCarthv fc et aJj |2008l. Unfortunately, 
these are the only reliable absolute fluxes available for our 
spectra. 

We observed J2123-0050 at the MMT observatory using 
the MMT Spectrograph in 2006.71 and 2009.40. The reso- 
lutions and wavelength coverages listed in Table 1 were ob- 
tained using a 1.5" wide entrance slit together with the 800 
groove/mm grating in the blue channel. We extracted and 
calibrated these spectra with standard techniques using the 
Low Redux Pipelin^] software package. We performed rel- 
ative flux calibrations using spectrophotometric standards 
observed on the same night. 

The emission-line redshift of J2123-0050 reported by the 



SDSS, 



2.2614, was m odified recently to 2.2686 ±0.0003 



by iHewett fc Wild! (|2010n . These redshifts match the peak 
position of the broad Civ AA1548,1551 emission line quite 
well. However, the lower ionization lines O I A1303 and C 11 
A1336 appear at a slightly higher redshift that is proba- 
bly closer to the true quasar systemi c (within ~100 km s _1 ; 
iTvtler fc Fanlll992l ; ls"hen et al'j|2007h . We therefore estimate 
the redshift from the O I and C n lines as follows. First we av- 
erage together the SDSS and MMT 2006.71 spectra to maxi- 
mize the signal-to-noise ratio. Then we manually interpolate 
across the tops of any narrow absorption lines to eliminate 
these features (e.g., in the Ol emission line profile), and cal- 
culate the 1 and C 11 emission line centroids. This yields an 
average redshift for the two lines of z e = 2.278±0.002 (using 



Year 


Date 


Telesc. 


Resol. 


AA 




(yyyy/mm/dd) 




(km s — 1 ) 


(A) 


2002.68 


2002/09/06 


SDSS 


150 


3805-9205 


2006.64 


2006/08/20 


Keck 


2.7 


3057-5896 


2006.71 


2006/09/16 


MMT 


220 


3050-4988 


2008.61 


2008/08/10 


VLT 


5.5 


3047-3870 




It 




4.5 


4622-9466 


2008.65 


2008/08/27 


VLT 






2008.66 


2008/08/30-31 


VLT 






2008.73 


2008/09/22 


VLT 






2008.73 


2008/09/25 


VLT 






2009.40 


2009/05/25 


MMT 


220 


3329-5968 



laboratory wavelengths from IVerner et al"1ll9961 ). Measure- 
ments based on the line peaks or the centroids of only the 
upper halves of the profiles yield essentially the same result. 
Also, redshifts determined separately for the two lines agree 
within the measurement uncertainty quoted above. We will 
adopt the redshift z e = 2.278±0.002 for the quasar systemic 
throughout the remainder of this paper. 

We observed J2123-0050 at high spectral resolutions 
at the W. M. Keck observatory using the High Resolu- 
tion Echelle Spectrograph (HIRES) and at the Very Large 
Telescope (VLT) using the Ultraviolet and Visual Echelle 
Spectrograph (UVES). We extracted and calibrated the 
Keck spectra using the HIRES reduction package HIRedu>Q. 
More information about the Keck /HIRES observati o ns an d 
data reductions can be found in IMilutinovic et all (|2010h . 
We obtained the VLT/UVES spectra in separate blue and 
red channels covering the wavelengths indicated in Table 
1. These data were reduced by the UVES data reduction 
pipeline. The final wavelength scales in both the Keck and 
VLT spectra are in the vacuum heliocentric frame. We com- 
bined the UVES spectra obtained in different echelle or- 
ders and in different exposures using the customized code 
UVES_poplei0. The VLT/UVES spectra were obtained on 
several dates during a ~6 week period in 2008 (Table 1). 
Careful examination of these spectra shows that there are no 
significant differences in any of the absorption lines within 
this data set. Hereafter, we discuss only the average VLT 
spectrum referred to by the average date, 2008.67. This spec- 
trum represents a total of 11.3 hrs of total exposure time on 
the quasar. 

To measure the absorption lines, we define pseudo- 
continua in all of the spectra by manually drawing a smooth 
curve that matches the data in every important detail while 
extrapolating over the tops of the absorption lines. We then 
divide the original data by these pseudo-continua to obtain 
normalized spectra for our analysis below. This procedure 
works well red ward of the the Lya emission line (A > 3985 
A observed) in all of our data. However, in the lower resolu- 
tion SDSS and MMT data, the blending of lines in the Lya 
forest is too severe to yield useful results. We therefore ex- 
clude the SDSS and MMT data at these short wavelengths 
from our analysis. 



1 http: //www. ucolick.org/~xavier/LowRedux/index. html 



2 http: / /www. ucolick.org/~xavier/HIRedux/index. html 

3 http://astronomy.swin.edu.au/~mmurphy/UVES_popler.html 
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2.1 Quasar Properties 

The photometric magnitudes reported by the SDSS (e.g., 
r = 16.4 and i = 16.3) indicate that J2123-0050 is one 
of the most luminous quasars known. (Inspection of the 
SDSS images reveals no evidence for gravitatio nal lensing 
that might artificially boost the observed fluxes. Ijust et all 
120071 ). We estimate the bolometric luminosity of J2123-0050 
to be L = 8.4 x 10 47 ergs s -1 based on the specific luminosity 
at 1450 A in the rest frame of AL A (1450A) « 1.9 x 10 47 ergs 
s _1 (see Figure 2 below) and a bolometric co rrection factor 
of L = 4.4 ALa(1450A) (see Appendix, also IWarner et al.1 
I2004T ). Combining this luminosity with our own measure- 
ment of the C IV emission line width, FWHM « 8750 km s _1 
from the SDSS spectrum, yie lds an estimate of the black hol e 
mass, Mbh ~ 2 x 10 10 M Q (|Vestergaard fc Petersor]|2006h . 
and the Eddington ratio, L/Le ~ 0.4, where Le is the Ed- 
dington luminosity. These values of Mbh and L/Le are 
within the nominal range of other high- luminosity quasars 
jShemmer et alJ^OO^lNetzer et al.ll2007fi . 

Some of our analysis in §3 and §4 below relies on size 
estimates for the quasar emission sources. If the UV con- 
tinuum source is a geometrically thin accretion disk that 
emits like a blackbody at every radius, then the size of the 
continuum emitter at different wavelengths depends only on 
the bolometric luminosi ty and the rad ial run of tempera- 
ture, T(r), in the disk |Petersonlll997h . The luminosit y of 
J2123-0050 combined with T(r) oc r -0 57 (|Gaskel]||200Sl ) in- 
dicates continuum source di ameters of ~0.026 pc at 1550 
A and ~0.013 pc at 1034 A (|Hamann fc Simorj|201ol ). The 
radius of the C IV broad emission line region (BELR) should 
be roughly Rbelr ~ 0.65 pc based on the e mpirical scaling 
rel ation with luminos i ty dBentz et al.l [2007) as formulated 
bv lHamann fc Simon! |2010l ). 

2.2 Other Data 

lJust et al.1 |2007h note that J2123-0050 is radi o quiet based 
on a non-detection in the FIRST radio survey (|Becker et al.1 
1 19951 ). They also report that J2123-0050 has a relatively 
"soft" X-ray spectrum, with a powerlaw slope of a < —1.1 
(/„ oc f a ) from roughly 0.5 to 8 keV (observed frame) and 
no evidence for continuous X-ray absorption in the quasar 
environment. The short exposure time of the X-ray obser- 
vations does not provide a strong upper limit on the col- 
umn density of X-ray absorption. However, the X-ray char- 
acteristics of J2123- 0050 are like oth er luminous radio-quiet 
quasars, for which lJust et al.1 (|2007l ) derive an upper limit 
of Nh < 2 x 10 21 cm~ 2 on the X-ray absorbing column 
using coadded spectra of quasars in their sample. Consis- 
tent with a lack of X-ray absorption, J2123-0050 has a 2- 
point power law index of a ox = 1.91 between 2500 A and 
2 keV in the quasar frame, which is in the nominal range 
for luminous radio-quiet quasars wi th no X-ray absorption 
jjust et al.ll2007l ; ISteffen et alj|2006h . 



3 ANALYSIS & RESULTS 

3.1 Defining the Variable Systems 

Figure 1 shows the five variable Civ AA1548,1551 doublets 
(labeled A-E) in the normalized Keck 2006.64 and VLT 



2008.67 spectra. Significant changes in the line strengths are 
clearly evident. It is also evident that the changes were well 
coordinated between the five systems (see also §3.2 below). 

Figure 2 shows the variable C IV NALs in relation to 
other features in the J2123-0050 spectrum. The spectrum 
shown in this figure represents the product of the normal- 
ized Keck 2006.64 data times the pseudo-continuum drawn 
through the SDSS 2002.68 spectrum. A small featureless gap 
in the Keck wavelength coverage around 4970 A (observed) 
is filled using the VLT spectrum. The result in Figure 2 thus 
has the fluxes, broad emission lines and overall shape match- 
ing the SDSS data in 2002.68, but with the absorption lines 
measured at much higher resolution at Keck in 2006.64. 

There are a number of other narrow C IV systems scat- 
tered across the J2123-0050 spectrum. They include one sys- 
tem at a smaller velocity shift than the variable lines (at 
~4960 A observed) and another that is blended with the 
variable system E (see also Figure 1). Careful comparisons 
between the Keck and VLT spectra shows that none of these 
other narrow C IV systems varied significantly (e.g., by more 
than a few percent in their rest equivalent widths) between 
the 2006.64 and 2008.67 observing epochs. Several Civ sys- 
tems measured in the Keck spectrum appear at wavelengths 
not covered by the VLT data. For these systems, we searched 
for variability in all of the other data sets listed in Table 1. 
We again find no evidence for variability in the other lines, in 
this case above an estimated sensitivity threshold of ~20% 
in the integrated line strengths. The only variable Civ lines 
detected in our data are the five systems labeled A-E in 
Figures 1 and 2. 

One important property of the variable C IV lines is that 
they all reach signficantly below the continuum shown by the 
dashed curve in Figure 2. We estimate this continuum based 
on a single powerlaw fit constrained by the measured flux 
in narrow wavelength intervals (1437-1448 A, 1675-1700 A 
and 1990-2030 A rest) that avoid the broad emission lines. 
The powerlaw index of the fit is a — —0.67 (for /„ oc 
The fact that the variable systems signficantly absorb the 
quasar continuum is important because the partial line-of- 
sight covering deduced from the line ratios (§3.3 and §3.4 
below) applies to the continuum source and not (or not only) 
to the much larger BELR (§2.1). 

Finally, we search for lines of other ions in the five 
variable systems using both the Keck 2006.64 and VLT 
2008.67 data. This search is complicated for wavelengths 
below ~1270 A in the absorber frame because of contam- 
ination by the dense forest of unrelated (intervening) Lya 
absorption lines. Nonetheless, we find strong absorption in 
Ovi AA1032,1038 in all five variable systems, probable ab- 
sorption in N v AA1239,1242 in some of those systems, and a 
secure detection of Lya in system C. We note that the detec- 
tion of O VI in system B confirms the reality of this system, 
which is not obvious from C IV alone because of the blend- 
ing with systems A and C (see Fig. 2 and Fig. 4 below). No 
other lines are detected in the variable systems, including 
specifically Ly/3 and low- ionization metal lines such as Cn 
A1336, Cm A977, Sim A1206 and Siiv AA1394,1403. 

3.2 Line Measurements &: Variability Properties 

Table 2 lists several parameters of the variable C IV lines 
A-E measured from the normalized Keck 2006.64 and VLT 
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Figure 1. Normalized spectra of J2123-0050 showing the variable Civ systems A-E in the two high-resolution observations: 2006.64 
(black curve) and 2008.67 (red). The velocity scale across the top applies to the short wavelength lines in the doublet, Civ A1548, 
relative to the emission redshift, z e = 2.278. The la variance spectra are plotted as dotted black and red curves across the bottom. Much 
narrower absorption lines not labeled are cosmologically intervening and unrelated to the variable outflow systems. 
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Figure 2. Synthesized spectrum of J2123-0050 showing the absorption lines in the Keck 2006.64 data relative to the broad emission 
lines and overall spectral shape defined by the SDSS measurement in 2002.68. All of the C IV AA1548,1551 absorption doublets detected 
in this wavelength range are marked by open brackets above the spectrum. The variable outflow systems are labeled A-E. Various broad 
emission lines are labeled across the top. A small featureless segment near 4970 A observed (drawn in red) uses the average VLT 2008.67 
spectrum to fill a small gap in the Keck wavelength coverage. The Keck and VLT spectra are shown after binomial smoothing to improve 
the presentation. The dashed curve is a powerlaw fit to the underlying quasar continuum. 
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2008.67 spectra. These results were obtained using cursor 
commands in the software package. For comparison, 

Table 2 also lists measurements (from the Keck spectrum 
only) for the three non- variable Civ systems that are clos- 
est in velocity to systems A-E. Column 1 in the table gives 
the variable system name, where A+B and C+B refer to the 
blended C IV features. Column 2 identifies the line within the 
C iv doublet. For the blends A+B and C+B, the doublet line 
listed pertains to systems A and C, respectively, which ap- 
pear to dominate the absorption based on our line fits in 
§3.3 below. Columns 3 and 4 give the observed wavelengths 
(Xobs, vacuum heliocentric) and redshifts (z a ) of the line cen- 
troids. Column 5 lists the centroid velocities (v) relative to 
z e = 2.278. Columns 6 and 7 give the rest equivalent widths 
(REW) and the full widths at half minimum (FWHM). The 
non-variable Civ systems at z a = 2.1291 and z a = 2.1139 
each consist of two or more blended sub-components. The 
table lists the total REWs for these blends while the FWHM 
pertains to the strongest single line in the blend. 

For the variable systems A and E, we record the REWs 
and FWHMs in Table 2 after removing contributions from 
unrelated non- variable lines (see Figs. 2, 4 and 5 below). 
This was achieved by extrapolating the A and E line profiles 
across the tops of the much narrower unrelated features. We 
are guided in this by the fact that the unrelated features did 
not vary between 2006.64 and 2008.67. Thus we can accu- 
rately decompose the blends into variable and non-variable 
components. 

In the lower resolution SDSS and MMT spectra, the 
C IV systems A+B+C and D+E are severely blended to- 
gether. We therefore measure only the total REWs in these 
blends for comparison to the Keck and VLT measurements. 
We also correct the total SDSS and MMT REWs for small 
contributions from the non-variable lines that are blended 
with systems A and E (as measured from the Keck and 
VLT spectra). Figure 3 plots the total C IV REWs in systems 
A+B+C and D+E at each observed epoch. The vertical bars 
in the figure are estimates of the la uncertainties associated 
with each REW measurement. These errors are dominated 
by the pseudo-continuum placement. We estimate the un- 
certainties by making repeated measurements of the REWs 
using different plausible continuum heights. 

One surprising result in Figure 3 is that the REWs 
differ dramatically between the Keck 2006.62 and MMT 
2006.71 spectra obtained only ~1 month apart in the ob- 
served frame. These REW differences are also evident from 
direct comparisons of the normalized spectra (not shown). 
They have a formal significance of 3<r for systems A+B+C 
and 6a for D+E (based on the uncertainties shown in Figure 
3) . We carefully reviewed our pseudo-continuum placements 
across these lines and find no evidence for uncertainties 
larger than those depicted in Figure 3. Nonetheless, we con- 
sider this evidence for dramatic short-term REW changes to 
be tentative because of the inherent difficulties in compar- 
ing the strengths of weak absorption lines in spectra with 
vastly different resolutions. In this case, the much shallower 



4 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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Figure 3. The total rest equivalent widths (REW) of the variable 
C IV systems A+B+C (filled red circles) and D+E (open black 
squares) are shown for each year observed. The vertical bars in- 
dicate approximate la uncertainties in the REW measurements. 



appearance of the blended lines in the lower resolution data 
makes them more susceptible to uncertainties in the pseudo- 
continuum placement that might be caused, for example, by 
small changes in the underlying emission line spectrum. We 
also note that significant line variations did not occur be- 
tween our much more precise VLT measurements obtained 
over a similar period of ~1.5 months in 2008 (Table 1). 

Figures 4 and 5 compare the normalized high-resolution 
Keck and VLT spectra across the variable C IV and O VI sys- 
tems. The O VI lines are severely blended with unrelated fea- 
tures in the Lya forest. Nonetheless, the Ovi outflow lines 
are clearly identified by their variability and their distinc- 
tive profiles and redshift agreement with C IV. Notice that 
the O VI lines are both stronger and broader than the C IV 
features, with variable absorption that appears to extend 
across a wider range of velocities compared to C IV. 

Figure 6 shows the Keck and VLT spectra across Lya 
in the variable systems A, B and C. Lya is clearly detected 
only in system C, based on its good match to the redshift, 
absorption profile and variability characteristics of C IV. Lya 
is significantly not present in system B, while its strength in 
the other systems (including D and E, not shown) is poorly 
constrained because of line blending problems. 

3.3 Covering Factors & Column Densities 

The smooth green and blue curves in Figures 4-6 are fits 
to the variable line profiles. These fits assume that at each 
velocity, v, a spatially homogeneous absorbing medium cov- 
ers a fraction, < C v ^ 1, of a spatially uniform emission 
source. In this situation, the line intensities seen by a distant 
observer are given by 

£ = (1-C v ) + C v e- T " (1) 

-ic 

where I c is the unabsorbed source intensity and r v and C v 
are the lin e optical depth and covering fraction, respec- 
tively fsee iGangulv et ail Il999l : IHamann fc Ferlandl 1 19991 ; 
IHamann fc Sabraf 2004 . for more general formulations). We 
can solve Equation 1 for t v and C v at each v using the 
intensity ratios of lines within multiplets, such as C IV 
AA1548,1551, where the optical depth ratio is known from 
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Table 2. Civ Line Measurements 



System 


Line 








2006.64 Keck 

V 

(km s J 










2008.67 VLT - 






^obs 

(A) 




Za 


REW 

(A) 


TIT. T TT TTV J 

FWHM 
(km s J 


^obs 

(A) 


Za V 

(km s- 1 ) 


REW 

(A) 


FWHM 
(km s ) 


A 


1551 


4921.3 


2 


17343 


-9716 





150 


62 


4921.3 


2.17342 


-9717 


0.063 


55 


A+B 


1548 


4913.2 


2 


17349 


-9710 





156 


66 


4913.2 


2.17348 


-9711 


0.061 


50 


C+B 


1551 


4905.0 


2 


16291 


-10710 





144 


85 


4905.0 


2.16295 


-10707 


0.055 


36 


C 


1548 


4896.7 


2 


16284 


-10717 





162 


62 


4896.7 


2.16285 


-10716 


0.082 


53 


D 


1551 


4865.0 


2 


13712 


-13162 





115 


100 


4865.3 


2.13731 


-13144 


0.024 


62 


D 


1548 


4856.5 


2 


13689 


-13184 





209 


164 


4856.6 


2.13694 


-13179 


0.075 


111 


E 


1551 


4850.6 


2 


12788 


-14045 





180 


147 


4850.4 


2.12775 


-14056 


0.114 


142 


E 


1548 


4842.4 


2 


12776 


-14056 





245 


154 


4842.8 


2.12799 


-14034 


0.133 


121 




1551 


4966.3 


2 


20246 


-6987 





024 


27.1 














1548 


4958.1 


2 


20247 


-6986 





042 


27.2 














1551 


4852.5 


2 


12909 


-13930 





077 


14.2 














1548 


4844.5 


2 


12910 


-13928 





123 


17.3 














1551 


4829.0 


2 


11395 


-15380 





090 


23.6 














1548 


4821.0 


2 


11394 


-15382 





163 


26.1 













the atomic physics, in this case T1548/T1551 ~ 2.0 (see also 
iHamann et al.ll 19971 ; fear low fc Sargenull997T ). For the reso- 
nance transitions discussed here, and assuming the ions that 
are entirely in their ground states, the ionic column densities 
are given by 



A^ ion 



Tve 2 fX 



dv 



(2) 



where / and A are the line o scillator strength and lab ora- 
tory wavelength, respectively (|Savage fc Sembacr][l99lh . 

Previous studies have shown that t v and C v can both 
have complex velocity-dependent behavio rs that differ 
signifi c antly between lines a nd bet we en ions (IHamann et all 
19971; iBarlow fe Sargentl Il997l; iGangulv et alj Il999l; 
Hamann et all l200ll; IHamann fc Sabral |2004| ; iGabel et al l 
2005l . l2006j ; lArav et al.ll2005l . 120081 ). Unfortunately, blending 
problems prevent us from obtaining D-dependent solutions 
to Equations 1 and 2 for most lines in the variable systems 
A-E. We therefore adopt a parameterized approach that 
involves i) fitting the C IV profiles with constant values of 
the covering fractions, C v — C , across the line profiles, and 
then ii) shifting and scaling the C IV fits to place constraints 
on other lines. This procedure also has the advantage 
of avoiding potentially large uncertainties in C v , t v and 
therefore Ni 0n that can be introduced by noise, e.g., in 
the line wings, when solving Equation 1 point-by-point at 
each v across the line profiles. The reliability of our fitting 
pro cedure is discussed furthe r at the end of this section (see 
also ISimon fc Hamannll2010i ) . 

Our fits assume gaussian t v profiles with potentially 
different values of the doppler width parameter on either 
side of line center. This functional form provides flexibility 
to account for some of the asymmetry in the measured lines, 
while still preserving a simple relationship between Ni on and 
the line center optical depth, r Q , namely, 



A^ ion 



3.34 x 10 1 



(3) 



where A D is the lab wavelength in A, and bb and b r are the 
blue and red side doppler parameters, respectively, in km 
s _1 . Each doublet pair is forced to have the same doppler 



parameters, velocity shift and a 2:1 ratio in the line optical 
depths. 

Our initial attempts to fit the C IV profiles using a x 2 
minimization technique led to spurious results. We found 
that manual adjustments and iteration by trial and error 
provide more flexibility and a much better understanding of 
the uncertainties in situations like this where the line profiles 
are irregular and the blending is potentially severe (e.g., in 
the Civ systems A+B+C). 

The fit results in Figures 4 and 5 show that the asym- 
metric gaussians provide a good match to the C IV profiles in 
systems D and E, but they are not able to produce simulta- 
neously the narrow cores and extended wings in systems A, 
B and C. We choose the "best" fits by giving highest priority 
to matching the line cores because they are most important 
for the accuracy of key parameters, r , C and Ni on . Note 
that the amount of C IV absorption contributed by the fully 
blended system B is highly uncertain. The line ratios across 
the blend A+B+C indicate that the amount of C IV absorp- 
tion in system B is less than A and C. Repeated attempts 
to fit C IV in the three systems A, B and C simultaneousljjf] 
indicate that system B might contribute up to ~30% of the 
total REW in the lines at this location. However, the data 
are also consistent with no C IV contribution from system B 
at all. Therefore none is included in the final fits shown in 
Figure 4 (also Table 3 below). 

Table 3 lists the ClV fit parameters derived from the 
Keck (2006.64) and VLT (2008.67) spectra. Columns 1 and 
2 indicate the system names and the observation years. 



5 There is no straightforward way to combine the intensities of 
blended lines with partial covering because the absorbing regions 
can have different values of C v across different spatial locations 
in front of the emission source. Our treatment assumes the maxi- 
mum amount of spatial overlap between absorbing regions. For 
example, if absorbing regions 1 and 2 have covering fractions 
Ci > Ci and line optical depths n and T2, respectively, then 
the total covering fraction is C\. The total line optical depth is 
ti + T2 in the spatial overlap region with covering fraction C2, 
and ti in the non-overlap r egion with covering fraction C\ — C2 ■ 
See lHall et al. I ||2003| , |2007|) for more examples and illustrations. 
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Figure 4. Normalized spectra showing the Civ (top panel) and Ovi (bottom) doublets in the variable systems A, B and C as observed 
in the Keck 2006.64 (black curve) and VLT 2008.67 (red) data. The doublet positions are shown by the brackets above the spectra. Many 
unrelated absorption lines are p resent in the O VI p lot caused by the Lya forest and by a system of extremely narrow H2 features in an 
intervening galaxy at z = 2.059 ( Malec et al. 2 01dl) . The velocity scale across the top of each panel applies to the short wavelength lines 
in the doublets, Civ 1548 and Ovi 1032, relative to the emission redshift, z e = 2.278. The smooth blue and green curves are fits to the 
variable lines in the Keck and VLT data, respectively, as described in §3.3. The dotted curves across the bottom are the ler uncertainties. 



Columns 3 and 4 give the redshifts, z , and velocity shifts, 
v , of the gaussian centers. (Note that these are not cen- 
troid positions unless the profile is symmetric.) Columns 5 
and 6 list the blue and red side doppler parameters. Column 
7 gives the line-center optical depths for the stronger A1548 
transition. Column 8 lists the covering fractions, and Col- 
umn 9 provides the ionic column densities from Equation 3. 
We estimate la uncertainties in these quantities by trial and 
error experiments with different fit parameters and different 
pseudo-continuum levels. The uncertainties in v , b r and bt 
should be ~3 km s _1 in systems A and C and roughly dou- 
ble that in the broader systems D and E. The covering frac- 



tions are in all cases well constrained, with errors of <10%. 
The optical depths and column densities have nominal la 
uncertainties around 25%. Including a maximum possible 
contribution from system B would lead to lower values of r 
and Ni 0n in systems A and C by about 30%. 

Next we transfer the C IV fits to other lines in the vari- 
able systems by keeping the kinematic parameters v , bf, and 
b r fixed while scaling only r and C , as needed. For system 
B, we adopt v — —10, 210 km s _1 and b r — bb = 30 km s _1 
based on our efforts to add a system B contribution to the 
C IV blend. This simple shift-and-scale approach is essential 
for dealing with non-detections and blends in the Lycv for- 
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Figure 5. Normalized spectra showing the Civ (top panel) and Ovi (bottom) doublets in the variable systems D and E. See Figure 4 
caption. 



Table 3. Civ Fit Results 



System 


Year 


Zo 


Vo 

(km s- 1 ) 


(km s _1 ) 


b r 

(km s _1 ) 


To 


Co 


N ion 

(10 13 cm" 2 ) 


A 


2006.64 


2.17350 


-9709 


30 


30 


2.3 


0.38 


15.7 


A 


2008.67 


2.17344 


-9715 


30 


25 


2.8 


0.19 


17.5 


C 


2006.64 


2.16278 


-10723 


30 


30 


2.7 


0.37 


18.4 


C 


2008.67 


2.16284 


-10717 


30 


35 


1.0 


0.20 


12.5 


D 


2006.64 


2.13727 


-13148 


95 


45 


1.5 


0.30 


23.8 


D 


2008.67 


2.13733 


-13142 


65 


10 


0.8 


0.18 


9.54 


E 


2006.64 


2.12790 


-14043 


80 


60 


2.8 


0.30 


44.5 


E 


2008.67 


2.12816 


-14018 


90 


50 


1.7 


0.22 


27.0 
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Figure 6. Normalized spectra obtained at Keck (2006.64, black 
curve) and the VLT (2008.67, red) showing the Lya line in the 
variable systems A, B and C. Other absorption lines, not labeled, 
are unrelated features formed in intervening clouds or galaxies. 
See Figure 4 caption. 



est. It is also justified by the good results shown for O VI and 
Lya in Figures 4-6. However, as noted in §3.2 above, the O VI 
absorption appears to span a wider range of velocities than 
C IV. Our fits ignore this additional O VI absorption because 
i) it is poorly measured due to blending problems, and ii) we 
want to limit our comparisons between C IV and O VI to the 
same kinematic gas components anyway, for the ionization 
and abundances analysis in §3.6 and §3.7 below. 

Table 4 lists the fit parameters derived for the other 
lines. The O VI doublets measured in 2006.64 and 2008.67 
provide useful constraints on both C and N . Notice that 
the Ovi covering fractions are larger than Civ (c.f., Table 
3). Also note that several of the O VI doublets are saturated, 
with intensity ratios close to unity (see Figs. 4 and 5). For 
these systems, Table 4 lists lower limits on the column den- 
sities that correspond to r = 4 in the stronger Ovi A 1032 
line. The values of Ni on listed for O VI system C and system 
E in 2008.67 (only) come from fits indicating r < 4. How- 
ever, we could conservatively consider all of the O VI column 
densities to be lower limits because of the blending problems 
in the Lya forest. We estimate lcr uncertainties in the Ovi 
results by repeating the fits with different continuum place- 
ments while keeping z and the doppler parameters fixed, 
as described above. We find that the errors in C are <15% 
while the errors in Ni 0n are roughly 30%. We conclude that 
the differences in the derived covering fractions between O vi 
and C IV, and the smaller changes in the O vi covering frac- 
tions compared to the changes in Civ (c.f., Tables 3 and 4), 
both have a large significance. 

The remaining lines listed in Table 4 are either non- 
detections (Si iv and C m) or severe blends (N v and H i) 
where we cannot use line intensity ratios to determine both 
To and Co- We therefore fix C along with the kinematic 
parameters and scale only r D to derive constraints on Ni 0n . 
Figure 6 shows examples of these fits applied to Lya in sys- 
tems A, B and C. We adopt values of C from C iv (Table 3) 
in all cases except for Lya in system B, which has no C IV fit 
result, and Nv, which has measured lines/blends too deep 



Table 4. Addtional Fit Results 



Line System Year Co Wj„ 

(10 13 cm" 2 ) 

Ovi 
Ovi 
Ovi 
Ovi 
Ovi 
Ovi 
Ovi 
Ovi 
Ovi 
Ovi 

Si iv 
Si iv 
Si iv 
Si iv 
Si iv 

Nv 
Nv 
Nv 
Nv 
Nv 

Hi 
111 
111 
Hi 
Hi 
Hi 

Cm 
Cm 
Cm 
Cm 
Cm 



to be consistent with C in C iv. For these latter cases we 
adopt the covering fractions from Ovi. Two of the entries 
for N v in Table 4 are blank because the blending is too se- 
vere to yield meaningful results. The remaining three Nv 
systems (B, C and E) are listed as upper limits in Table 4 
because of blending uncertainties. We note, however, that 
there do appear to be real Nv lines present with redshifts 
and profiles similar to C iv. The upper limits on iVj on listed 
for Hi in Table 4 use the best available constraints from 
either Lya or Ly/3. The specific values of the Hi column 
density given for system C come from the Lya fits shown 
in Figure 6 without additional constraints from Ly/3. The 
errors in these H I results might be as large as a factor of ~2 
due mainly to the uncertainties in C . 

As a further check of our fitting analysis, we derive 
velocity-dependent values of r v and C v for the two best- 
measured cases: the C iv lines in systems D and E observed 
in 2006.64. Specifically, we calculate the average line inten- 
sities in velocity bins 25 km s _1 wide and then solve Equa- 
tion 1 for t v and C v in each bin. Figure 7 compares the 
derived values of 1 — C v to the observed line profiles. The 
vertical bars running through the filled circles in Figure 7 
represent la uncertainties in C v caused by photon statis- 



A 2006.64 

A 2008.67 

B 2006.64 

B 2008.67 

C 2006.64 

C 2008.67 

D 2006.64 

D 2008.67 

E 2006.64 

E 2008.67 

A 2006.64 

B 2006.64 

C 2006.64 

D 2006.64 

E 2006.64 

A 2006.64 

B 2006.64 

C 2006.64 

D 2006.64 

E 2006.64 

A 2006.64 

B 2006.64 

C 2006.64 

C 2008.67 

D 2006.64 

E 2006.64 

A 2006.64 

B 2006.64 

C 2006.64 

D 2006.64 

E 2006.64 



0.95 >58.4 

0.82 >53.5 

0.38 >58.4 

0.30 >58.4 

0.87 51.1 

0.82 55.4 

0.50 >136 

0.50 >102 

0.54 >136 

0.63 92.0 

<0.30 
<0.30 
<0.30 
<0.60 
<0.60 



<8.5 
<18 

<72 

<5.2 
<0.6 
4.8 
7.7 
<6.4 
<23 

<5.0 
<3.0 
<3.0 
<15 
<15 
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Figure 7. C IV line profiles measured in 2006.64 in systems D (top 
panel) and E (bottom), compared to derived values of 1— C v (filled 
circles), on a velocity scale relative to z e = 2.278. The solid and 
dashed curves show the stronger (A1548) and weaker (A1551) dou- 
blet members, respectively, after smoothing the spectra 3 times 
with a binomial function. The vertical bars on the filled circles 
are la uncertainties. For the velocity bins with r v -C 1 (in the 
blue line wings), two filled circles at each velocity represent the 
full range of 1 — C' v values consistent with the data. 



tics (|Hall et al.ll2003h . However, in the blue wings of both 
systems D and E there are several velocity bins where the 
doublet ratios indicate r v <C 1 and the standard analysis 
yields C v = 1. The actual values of r„ and C v are poorly 
constrained in these cases because the right side of Equation 
1 reduces to 1 — C v t v in the t„<1 limit and, therefore, the 
line intensities constrain only the product C v t v instead of 
C v and t v separately. For example, in the velocity bin cen- 



tered at v 



-13,250 km s in Civ system D, the data 



are consistent with solutions ranging from 1 — C v = and 
t„(1548) = 0.14 to 1 - C v « 0.5 and t v (1548) = 0.28. Fig- 
ure 7 shows both the maximum and minimum (1 — C v =0) 
results, connected by a vertical bar, for each of the optically 
thin data points. The data at these velocities are compatible 
with any value of 1 — C v within the ranges shown. 

The values of C v and r„ derived from the point-by-point 
analysis near line center agree well with C and r obtained 
above from the gaussian fits (Table 3). The point-by-point 
results in Figure 7 show further that G v changes significantly 
across the line profiles. In particular, the shape of the inten- 
sity profile across the core of system E appears to be gov- 
erned largely by the behavior of C v with t v > 2. In system D, 



there appears to be an abrupt change in C v around —13,200 
km s _1 , between a line core that has t„ > 1 and C v « 0.28 
and a blue wing with t v <C 1 and C v ~ 1. Integrating t„ over 
all of the velocity bins shown in Figure 7 (as per Equation 
2 but excluding the lowest v point in the red wing of both 
systems because they formally yield t„ 1 due to noise) 
indicates Ni 0n = 2.3 ± 0.3 x 10 14 cm -2 for system D and 
3.8 ± 0.3 x 10 14 cm -2 for system E. These column densities 
are in excellent agreement with the gaussian fit results in 
Table 3. 

We conclude with a comment about the assumption of 
spatial uniformity in the absorbing regions. Some studies of 
partial covering in quasar absorption lines indicate that the 
the actual absorbing structures are inhomogeneous, with a 
range of r„ values at each v across the projected area of the 
emission source (lHamann et al.l l200ll; Ide Kool et alj |2002| ; 
lHamann fc Sabrall2004l ; lArav et al.ll2008l . and refs. therein) . 
The evidence for this is significantly different covering frac- 
tions in different lines/ions. While these data can provide 
important constraints on the nature of the absorber inho- 
mogeneities, there is not enough information to determine 
uniquely the 2-dimensional spatial distributions of t v in 
each line. The situation is further complicated by the fact 
that the background light sources, e.g., the accretion disk 
and/or the BELR, have their own spatial non-uniformities 
that can include a wavelength-dependent size. A good strat- 
egy, therefore, is to consider the observationa l consequences 
of different plausible t v spatial distributions ( de Kool et al.l 



|2002| ; lHamann fc Sabrall2004l ; lArav et al.ll20051 . l2008l ). In one 

study, extensive simulations using different functional forms 
for the t v spatial distribution show that Equation 1 (which 
assumes a top hat t„ distribution) yields good approxima- 
tions to the spatially-averaged values of t v and Nj on in more 
gener al inhomogeneous absorbing media |Hamann fc Sabral 
120041 ). Therefore, the results derived from Equation 1 are 
appropriate for making spatially-averaged estimates of the 
ionization and metal abundances (§3.6 below). They are also 
useful for characterizing the covering fractions even for com- 
plex absorbers. In particular, the values of C v derived from 
Equation 1 correspond roughly to the amount of coverage 
by absor bing material with t v > 1 in an inhomogeneous 
medium l|Hamann fc Sabralf2004h . 

3.4 Absorbing Region Sizes 

The partial covering discussed in §3.3 implies that the ab- 
sorbing structures are not much larger than the background 
light source. (The absorbers could be smaller than the emis- 
sion source, but partial covering by much larger clouds would 
require unphysically sharp edges and extremely unlikely spa- 
tial alignments.) Figure 2 shows that the variable C IV sys- 
tems lie in the far blue wing of the C IV BEL, where the flux 
is dominated by continuum emission. In principle, the line 
absorption could occur in front of both the UV continuum 
source and the BELR, with potentially diff erent covering 
fractions for each region jGangulv et alll999l) . However, the 
depths of the lines below the quasar continuum (Figure 2) 
imply that the C IV partial covering results pertain mostly or 
exclusively to the continuum source. This conclusion is con- 
firmed by the Ovi lines, which have even deeper absorption 
troughs and no significant underlying BELR flux. 

The extent to which the absorbers cover the BELR are 
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thus undetermined by our data. This ambiguity introduces 
a small additional uncertainty in C v derived for the C IV 
lines because the BELR contributes roughly 10% of the to- 
tal flux at those wavelengths (Figure 2). Our analysis in 
§3.3 assumed implicitly that C v is the same for all emission 
sources beneath the absorption lines. However, if the BELR 
is not covered at all, the continuum covering fractions would 
be roughly 10% larger than the values listed for C IV in Ta- 
ble 3. Conversely, if the C IV absorber completely covers the 
BELR, the continuum covering fractions would be ~20% 
smaller than the res ults listed in Table 3 (see equations in 
iGangulv et al.l 1 19991 ) . This latter situation (complete cov- 
erage of the BELR with partial covering of the continuum 
source) seems highly unlikely given that the C IV BELR is 
of order 20 times larger than the continuum source at 1550 
A (§2.1). Hereafter, we will assume that all of the derived 
covering fractions apply only to the continuum source. 

Comparing the covering fractions in C IV and O VI (Ta- 
bles 3 and 4) to the theoretical diameters of the continuum 
source at 1550 A and 1034 A, respectively (§2.1), indicates 
that the C IV and O VI absorbing regions both have charac- 
teristic sizes of order ~0. 01-0. 02 pc (as measured in 2006.64). 
It is important to keep in mind, however, that this charac- 
teristic size is actually just an upper limit on the true sizes of 
absorbing structures in the flow. In particular, partial cover- 
ing might be caused by patchy dist ributions of many smalle r 
sub-structures (see §4.4 below, also lHamann fc Sabrall2004l ). 

3.5 Outflow Dynamics 

The five distinct outflow systems span a range of veloci- 
ties from —9710 to —14,050 km s _1 and Civ line widths 
from FWHM « 62 to 164 km s _1 (Table 2). The Ovi lines, 
representing higher ionization material, appear somewhat 
broader than C IV in terms of both their FWHMs and broad 
weak wings (see Figs. 4 and 5). The outflow systems have a 
natural division into two groups, A+B+C and D+E, based 
on similarities in the line strengths, profiles and velocity 
shifts (see also Fig. 1). One remarkable characteristic of 
these systems overall is their small velocity dispersions com- 
pared to the high flow speeds, e.g., FWHM/?; ~ 100. This 
in marked cont rast to BAL outflows , which typically ha ve 
FWHM/w ~ 1 ijWevmann et al.lll99ll : Ikorista et al J 1 19931 ). 

The kinematic stability of the narrowest systems A, B 
and C between 2006.64 and 2008.67 places an upper limit 
on the radial acceleration of <3 km s _1 yr _1 in the quasar 
frame. We obtain this estimate from the measurements in 
Table 3 and by careful comparisons of simultaneous fits to 
all three narrow systems (in Civ and Ovi) in both epochs. 
(The centroids of the broader systems D and E did change 
slightly between observations, but this appears to be the re- 
sult of profile changes rather than a true velocity shift; Fig- 
ure 5.) For comparison, the gravitational acceleration due to 
Mbh at the radius of the Civ BELR (§2.1) is roughly 200 
km s" 1 yr _ . Therefore, the outward driving force is either 
in a delicate balance with gravity or, more likely, the ab- 
sorbing gas is located in a low-gravity environment farther 
from the black hole than the C IV BELR. A radial distance 
of R > 5 pc is needed for the local gravity to be less than the 
observed acceleration limit. At this distance, the observed 
outflow speeds vastly exceed the local gravitational escape 
speed. Thus it appears that the absorbing gas is gravitation- 



ally unbound and coasting freely at speeds near the outflow 
terminal velocity. 

There is an apparent occurrence of "line-locking" be- 
tween the narrow C IV systems A, B and C. In particular, 
the velocity differences between systems A-B and B-C are a 
close match to the C IV doublet separation of Av — 498 km 
s . (Although the system B lines are entirely blended in 
C IV, their velocity is confirmed by our fits to the O VI lines 
shown in Figure 4.) 

Line-locking is u sually interpreted as evidence for ra- 
diativ e acceleration |Foltz et all Il987 1_ Braun fc Milgroml 
I l989l: ISrianand fc Petitieanl l2000l : ISrianand et alJ \200% . 
IGangulv et al.l 12003') . Cases of multiple line-locks between 
systems with narrow profiles, as in J2123-0050, provide 
the best evidence for physically locked systems because 
they are difficult to explain by chance alignments of lines 
in the spectrum. However, it is not known if radiative 
forces can actually lock doublets like Civ together in real 
qu asar outflows. The most plausible scenario was outlined 
bv lBraun fc Milgroml (|l989l ): Consider two blobs of gas mov- 
ing along the same radial path in a radiatively accelerated 
outflow. The faster blob is farther from the emission source 
and experiencing a greater outward acceleration. When the 
velocity difference between the blobs matches the C IV dou- 
blet separation, Av m 498 km s~\ the A1551 line in the 
outer blob falls in the shadow cast by A1548 in the slower 
inner blob. This leads to diminished acceleration in the outer 
blob and a velocity lock between the two blobs at the dou- 
blet separation. Line-locking might even contribute to the 
formation of blobs because the localized dip in the radia- 
tive force downstream from the inner blob creates a n atural 
collecting area for outflow gas (| Ganguly et al. 2003). This 
scenario does not require that the C IV lines dominate the 
overall acceleration, only that the acceleration transferred 
to the outer blob via the A1551 line is at least as large as 
the acceleration difference between the two blobs in the ab- 
sence of shadowing (i.e., outside of the Av « 498 km s _1 line 
overlap situation). The appearance of a double line-lock in 
Civ, as in J2123-0050, might also be facilitated by the Nv 
absorption doublet whose separation, Av = 964 km s _1 , is 
close to the velocity difference between systems A and C. In 
any case, the feasibility of this line-lock interpretation has 
not been demonstrated in any quasar outflow. 

Here we note simply that if radiative forces did indeed 
lock the Civ doublets in J2123-0050, then the absorbing 
clouds must be moving almost directly toward us. This is 
because the velocity difference between line-locked blobs 
should match the C IV doublet separation along the flow 
direction. If the flow is not aimed directly at us, we should 
observe velocity differences between the systems that are 
smaller than the doublet separation. Our measurements of 
the line centroids (Table 2) and our fits to the profiles, in- 
cluding Ovi in system B (Table 3 and Figure 4), indicate 
that the velocity differences between systems A, B and C 
are within ~20 km s _1 of Av = 498 km s _1 . This implies 
that the absorber trajectories are within ~16° of the purely 
radial (line-of-sight) direction. 

3.6 Physical Conditions 

The strong detections of O VI compared to C IV absorption, 
and the non-detections of lower ionization lines like Si IV and 
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Cm (Tables 3 and 4), indicate that the degree of ionization 
is high and approximately the same in all of the variable sys- 
tems. The Appendix bel ow describes numer ical simulations 
using the code Cloudy (|Ferland et al.l fl998) of the ioniza- 
tion and physical conditions in clouds that are in photoion- 
ization equilibrium with a quasar continuum source. Com- 
paring those results (Figure Al) to the most stringent limit 
on the observed AT(ClIl)/iV(Clv) < 6 column density ratio 
(system C in 2006.64) indicates that the dimensionless ion- 
ization parameter is at least log U > —0.8. The measured 
ratios of N(0 Vl)/iV(C iv) > 3 provide a very similar con- 
straint on log U if we make the additional assumption that 
the O/C abundance is approximately solar. 

This lower limit on the ionization combined with the 
average measured Hi column density in system C, 7V(Hl) ~ 
6 x 10 13 cm -2 (Table 4), implies that the minimum total 
hydrogen column density in that system is Nh > 3 x 10 18 
cm~ 2 (see Figure Al). In §4.3 below, we argue that the line 
changes were caused by changes in the degree of ionization 
and the most likely value of the ionization parameter is 
log U ~ —0.4, near the peak in the O VI ion fraction. Us- 
ing that result, our best estimate of the actual total column 
density in system C is of order Nh ~ 10 19 cm -2 . Given the 
similarities between the systems A-E, we crudely estimate 
that the total column density in all five outflow absorbers is 
N H ~ 5 x 10 19 cm -2 . 

The more conservative requirement for log U > —0.8 
in a photoionized gas at R > 5 pc (§3.5) from the quasar 
places an upper limit on the volume density of nn < 2 x 10 8 
cm -3 (see Equation Al). A lower limit on the density can 
be derived by noting that, in photoionization equilibrium, 
changes in the ionization state require at least a recombi- 
nation time given by t r ~ (ariij)" 1 , where n e is the elec- 
tron density and a r is the recombination rate coefficient. 
Given a r « 10 -11 cm 3 s" 1 for either Ovn — s> Ovi or Ovi 
— > O v at a nominal gas temperature of T e ~ 20, 00 K 
fsee lArnaud fc Rothenfludll985l ; lHamann et al.lll995h , the 
observed Ovi variability in ^0.63 yr indicates a minimum 
density of n e ~ ng > 5000 cm~ 3 . This density combined 
with the ionization constraint log U > —0.8 yields a maxi- 
mum distance of 7? < 1.1 kpc (Eqn. Al). 

The density and distance constraints based on the re- 
combination time (riH > 5000 cm -3 and R < 1.1 kpc) might 
not apply if the ionization changes occurred out of equilib- 
rium. In particular, ionization increases caused by a rise in 
the ionizing flux could occur much faster than the recombi- 
nation time. The most likely explanation for the line changes 
between 2006.64 and 2008.67 is, indeed, an ionization in- 
crease (§4.3 below). However, there is weak evidence for both 
weakening and strengthening of the C IV lines in J2123-0050 
(Fig. 3), which is typical variable outflow lines in quasars 
lHamann et all Il995 ; IMisawa et al.ll2007l; iNaravanan et all 
20041; lHamann et al Il997l; iRodriguez Hidalgo et all l20ld ; 
Gibson et all l20ld ; ICapellupo et alJl201Ch . If this up and 
down behavior of the line strengths repr esents changes in 
the ionization (see also lMisawa et akll2007t ). then recombina- 
tion must be involved. Continued high-resolution monitoring 
is needed to confirm that conclusion specifically for J2123- 
0050. However, even if the line variations in this object oc- 
curred because of clouds crossing our view of the continuum 
source, basic assumptions about the crossing speeds would 
require a radial distance much smaller than the 1.1 kpc limit 



derived from the recombination time (§4.3). Therefore, we 
adopt R < 1.1 kpc as a firm upper limit. 

3.7 Metal Abundances 

We estimate the C/H metal abundance from the average 
measured ratio of iV(C Iv)/iV(H i) ~ 2.5 in system C (Tables 
3 and 4) and an ionization correction based on log U ~ —0.4 
(see Eqn. A2 and Fig. Al). This indicates [C/H] ~ 0.3 ±0.2, 
where the square brackets have their usual meaning of log 
abundance relative to solar and the la uncertainty includes 
an estimate of the errors in both the column densities and 
the ionization correction. The uncertainties related to the 
ionization are small because the most likely ionization pa- 
rameter, log U ~ —0.4, is near the peak in the C IV ion frac- 
tion curve where the correction factor, /(H i)//(C iv), is not 
sensitive to the specific value of log U across a fairly wide 
range. A careful inspection of the numerical results depicted 
in Figure Al shows that /(H i) //(C iv) stays constant within 
±0.1 dex across a range in logt/ from —1.4 to 0.0. More- 
over, this range in log U encompasses the minimum value 
of /(Hi)//(Civ) at \ogU = -0.75 (Table Al). If we apply 
the minimum correction factor log(/(H l)//(C iv)) = —3.80 
(Table Al) to the system C column densities, we derive a 
lower limit on the metal abundance, [C/H] > 0.2, which is 
close to our best estimate above but co mpletely indepe ndent 
of the ionization uncertainties (see also lHamannlH997l ) . 



4 DISCUSSION 

Section 3 presented measurements and analysis of five dis- 
tinct outflow NAL systems in J2123-0050 having Civ line 
widths of FWHM ~ 62 to 164 km s" 1 and velocity shifts 
from v ~ —9710 to —14,050 km s _1 in the quasar rest frame. 
These systems appear to be physically related based on their 
roughly similar line strengths, kinematics, ionizations, line- 
of-sight covering factors, and coordinated variabilities. They 
provide a wealth of information about a particular NAL out- 
flow that is valuable for comparison to other work on BAL 
and mini-BAL outflows. Here we discuss some additional 
results and implications. Please see §5 below for an overall 
summary. 



4.1 Location Summary 

The formation of the variable NALs in a quasar outflow is 
confirmed by all three properties of line variability, partial 
covering and resolved profiles that are smooth and broad 
compared to thermal line widths. Moreover, if we accept 
that the lines form somewhere in the quasar environment 
based on these properties, then the only plausible location is 
a quasar-driven outflow because the velocity shifts are much 
too large to be explained by other near-quasar environments 
such as a galactic starburst-driven w ind or a nearby galaxy 
in the same cluster as the quasar |Prochaska fc Hennawil 
120091 ; iRupke et al.l [2005V By matching the variability to a 
recombination time, we place an upper limit on the distance 
between the absorber and the quasar of R < 1.1 kpc (§3.6). 
At the opposite extreme, constraints on the acceleration 
suggest that the NAL gas is coasting freely at a distance 
R > 5 pc (§3.5). Thus the full range of plausible distances is 
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5 iS R iS 1100 pc. However, if the absorbing structures are 
small blobs or filaments created in the inner flow, their most 
likely location is near the ~5 pc minimum radius because 
such structures will travel just a few pc before dissipating 
(if there is no external confinement. §4.4 below). 

Finally, our estimate of the metallicity, [C/H] ~ 0.3 ± 
0.2 with a lower limit of [C/H] > 0.2, is also consis- 
tent with the formation of these NALs in a quasar out- 
flow. Super-solar metallicities are extremely rare in in- 
tervening absorption line systems, having been observed 
so far o nly in a few Lyman-limit and super-Lyman -limit 
systems i|Prochaska et al.l 120061 ; iProchter etafl feoiCh , but 
they are typical of near-quasar enviro nments based on es- 
timates from the broad emission lines (lHamann fc Ferlandl 
ll999l : lDietricn" et al ] |2003l ; lNagao" et al ]|200rj) and from other 
well-studied cases of NAL outflows (ID' Odorico et alj|2004 
iGabel et all 120051 . l200rj ; lArav et al l 12001 Simon. Hamann 

6 Pettini 2010, in prep.). 



4.2 The Outflow Origins of Quasar NALs 

It is interesting to note that the outflow NALs in J2123-0050 
are unresolved and indistinguishable from cosmologically in- 
tervening lines in medium resolution spectra like the SDSS. 
They are much too narrow to be identified with an outflow 
using on e of the absorption line indices designed for this 
purpose (|Wevmann et al.lll99ll ; iTrump et all 120061 ). Their 
velocity shifts are also well above the nominal \v\ < 5000 
km s _1 cutoff used to define "associated" absorption lines 
(AALs), which are likely to have a physical relationship to 
the quasar based o n their statistical exce s s at z a z e com- 
pared to z a < z e i Wevmann et al.lll979l ; iFoltz et al.l 1 19861 ; 
iNestor et all 120081 ; Iwild et al J 12008 ). Without high resolu- 
tion spectra or multi-epoch observations to test for variabil- 
ity, the outflow NALs in J2123-0050 would be mistakenly 
attributed to cosmologically intervening gas. 

The outflow systems are, in fact, surrounded in the 
J2123-0050 spectrum by other narrow C IV lines that did 
not vary, do not have partial covering and almost certainly 
do form in unrelated intervening material (Figure 1). One of 
these intervening systems is at a smaller velocity shift than 
the outflow lines (at v ~ 6987 km s _1 ), while another is 
blended directly with the outflow system E (see also Table 
2 and Figure 5) . This mixture of lines in the same spectrum 
clearly demonstrates that velocity shift alone is a poor in- 
dicator of the outflow v ersus intervening origin of narrow 
absorption lines (see also ISimon fc Hamannll20ld ). 

The secure outflow origin of the NALs in J2123-0050 
lends anecdotal support to recent claims that a significant 
fraction of narrow C IV systems, even at large velocity sh ifts, 
form in quasar outflows. For example, Richards (2001) ar- 
gued that ~36% of high- velocity NALs originate in outflows 
based on a correlation between the numbers of these systems 
detected and the rad io properties of the background quasars. 
iMisawa et~ai1 (|2007l ) estimated that 10-17% of narrow high- 
velocity C IV systems belong to quasar outflows based on 
the incidence of partial covering measured in high-resolution 
spectra. Similar results we re obtained more re cently by Si- 
mon et al. (2010, in prep.). lNestor et al.l (|2008l ) showed that 
the statistical excess of Civ NALs near the quasar redshift 
extends out to at least ~12,000 km s~\ e.g., well beyond 
the nominal cutoff of 5000 km s _1 for "associated" absorp- 



tion lines. They argue that >43% of NALs stronger than 
REW = 0.3 A in the velocity range 7 50 < v < 12, 000 km 
s _1 originate in quasar outflows (also Wild et al. 20081) . In 
their analysis of mostly weaker systems, Simon et al. (2010, 
in prep.) find that roughly 20-25% of C IV NALs in this ve- 
locity range exhibit partial covering. 

The variable lines in J2123-0050 are unique in that they 
are (so far) the highest velocity NALs known to originate in 
a quasar outflow based on all three indicators of line vari- 
ability, partial covering and super-thermal line widths. Their 
closest analogues in the literature are a complex of NALs at 
somewhat lower speeds, v ~ 9500 km s~\ in the re dshift 
2.5 quasar HS 1603+3820 (jMisawa et al.ll2005l , [20071 ) . and 
an isolated broader (FWHM ~ 400 km s _1 ) outflow system 
at v ~ 24,000 km s" 1 in the z e « 2.5 quasar Q2343+125 
l|Hamann et al.lll997l ). 

4.3 What Caused the NAL Variability? 

The key properties of the NAL variability in J2123-0050 
are 1) the line variations were well coordinated between 
the five systems, each with larger changes in C IV than 
Ovi, 2) changes in the NAL strengths were accompanied 
by roughly commensurate changes in the line-of-sight cover- 
ing fractions, and 3) the variability time is ^0.63 yr in the 
quasar rest frame. 

Changes in the covering fractions would seem to have 
a natural explanation in clouds moving across our lines 
of sight to the background emission source. This expla- 
nation is, i n fact, favored by s ome recent studies of BAL 
variability ilGibson et all 120081 . |2010| ; lHamann et~afl 120081 ; 
ICapellupo et alj2010|) . However, in J2123-0050, line changes 
induced by crossing clouds appear unlikely because the co- 
ordinated variations in five distinct NAL systems would 
require highly coordinated movements bet ween five dis- 
tinct absorbing structures in the outflow (also IMisawa et al.l 
2005). Moreover, crossing clouds would need to traverse a 
significant fraction of the emission source diameter in ^0.63 
yr to cause the observed line variations. In §3.4 we argued 
that the relevant emission source is the accretion disk with 
diameter -D1550 ~ 0.026 pc at 1550 A. Crossing ^15% of that 
length in ^0.63 yr, to explain the observed changes of order 

0. 15 in the Civ covering fractions, would require transverse 
speeds v tr > 6, 000 km s _1 . This seems unrealistic given that 
the radial flow speeds are not much larger, e.g., v ~ —9710 
to 14,050 km s" 1 . The maximum transverse speeds we might 
expect in a quasar outflow launched from a rotating accre- 
tion disk are of order the virial speed at the absorber loca- 
tion. A virial speed consistent with v tr > 6000 km s _1 in 
J2123-0050 would require that the absorber resides at a ra- 
dial distance <2.4 pc from the SMBH. This would place the 
absorbers inside of the minimum radius of 5 pc estimated 
above from the lack of radial acceleration (§3.5). It would 
also contradict the evidence from line-locking that the flow 
we observe is traveling primarily in a radial direction (§3.5). 

A more likely explanation for the coordinated line vari- 
ations is global changes in the ionization caused by fluc- 
tuations in the quasar's continuum flux. This situation 
could produce the observed changes in the covering frac- 
tions (§3.3) if the absorbers are spatially inhomogeneous, 

1. e., with a distribution of column densities across the emis- 
sion source, such that the projected areas with t v > 1 in 
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a given line change wi th the ionization state of the gas 
|Hamann fc Sabrall20o4) . The observed decrease in the C IV 
line strengths between 2006.64 and 2008.67, accompanied by 
smaller or negligible changes in O VI, are consistent with an 
overall increase in the ionization parameter if the starting 
point was near the peak in the O VI ion fraction. In par- 
ticular, an ionization parameter initially near log U ~ —0.4 
in 2006.64 and then increasing by a few tenths of a dex by 
2008.67 would cause a decline in the C IV line strengths and 
column densities consistent with the observations (Table 3), 
while the Ovi remains nearly constant (see Figure Al). 

This scenario requires that the quasar's far-UV con- 
tinuum flux changed by a factor of roughly two in ^0.63 
yr (rest). Little is known about the far-UV variability 
properties of luminous quasars. Nonetheless, this type of 
continuum change seems plausible given that 1) near-UV 
flux variations as large as 30-40% have been observed in 
other luminous q uasars on times scales of months (rest, 
iKaspi et alj|2007h . and 2) flux changes in the far-UV could 
be much larger than the near-UV because shorter continuum 
wavelengths generally exhibit larger amplitu de fluctuations 
jKrolik et al.lll99ll ; IVanden Berk et alj|20u4 ). 

The recent SDSS stripe 82 survey shows specifically 
that the near-UV continuu m variability of J21 23-0050 is 
typical of luminous quasars (|Schmidt et al.| [2010V This sur- 
vey includes roughly 60 epochs of imaging photometry col- 
lected over a ~7 year period that overlaps with our spectro- 
scopic observations. Examination of these data for J2123- 
0050 (kindly provided by Kasper Schmidt) shows contin- 
uum variations up to ~0.27 magnitudes (or ~28%) on ob- 
served time scales of ~2 years in the r band (~2000 A rest). 
Moreover, the r-band flux increased by this amount between 
roughly 2005.8 and the last stripe 82 data point in 2008.0, 
which is consistent with an our conclusion that the ioniza- 
tion increased between the Keck and VLT observations in 
2006.64 and 2008.67. 

The main caveat to interpreting the NAL variabilities 
in terms far-UV flux changes is that there might have been 
a large NAL variation between 2006.64 (Keck) and 2006.71 
(MMT), i.e., in just ~8 days in the quasar frame (Figure 
3). If this extreme short-term variability can be confirmed, 
it might be difficult to explain via factor of ~2 changes in 
the far-UV emissio n from the accretion disk (but see also 
iMisawa et"afll2007fi . 



4.4 Outflow Structure 

The outflow NALs in J2123-005 require at least five distinct 
absorbing structures that have similar physical conditions 
(ionization, column densities, covering fractions, §3.3 and 
§3.4) and roughly similar kinematics (lines dispersions and 
velocity shifts, §3.5). The partial covering evident in all five 
systems indicates that the absorbers are each characteristi- 
cally smaller than the quasar's UV continuum source, e.g., 
~0. 01-0. 02 pc across. All five absorbers must also be spa- 
tially inhomogeneous to explain the observed changes in the 
covering fractions if those changes are, indeed, caused by 
changes in the quasar's ionizing flux (§4.3). 

The absorbing structures might be discrete blobs or fil- 
aments created by instabilities in the inner flow, perhaps 
resembling the clumpy str uctures seen in some nu merical 
simulations of BAL winds jProga fe KaHmanll2004h . If this 



assessment is correct, then their most likely location is near 
the minimum radius of 7? ~ 5 pc inferred from the lack of 
acceleration (§3.5). This is because small blobs or filaments 
have short survival times against dissipation if there is no 
external pressure confinement. The characteristic absorber 
diameter of d ~ 0.015 pc (§3.4) combined with the doppler 
parameters b ~ 70 km s _1 of the broader systems D and E 
(Table 3) indicate dissipation times of order tdis ~ d/b ~ 210 
yr. This is much longer than our observation timeline (Table 
1) and therefore not a factor in the line variability. How- 
ever, gas blobs moving radially at v ~ 13, 500 km s _1 (for 
systems D and E) will travel a distance of just ~3 pc be- 
fore dissipating. Therefore, they should be located within 
~3 pc of their point of origin. If that was the inner flow 
near the accretion disk (e.g., at rad ii not much greater than 
a few times Rbel r ~ 0.65 pc, §2.1. lMurrav fc Chiand[l997l ; 
iProga et ai]|200Ch . then the absorbers should reside near the 
R ~ 5 pc minimum radius (§4.1). 

One puzzling constraint on the internal structure of 
these absorbers comes from the high densities. At R ~ 5 pc, 
the densities need to be near the upper end of the estimated 
range, 10 8 > uh > 5000 cm -3 , to avoid over- ionization by 
the intense quasar radiation field (§3.6). However, this entire 
density range is much larger than the average gas density, 
(riff) ~ 215 cm" 3 , expected for absorbing clouds with total 
column density Nh ~ 10 19 cm -2 (§3.6) and characteristic 
size ~0.015 pc. This disparity implies that each of the five 
absorbers is composed of dense sheets or sub-structures that 
have radial thickness much smaller than their characteristic 
transverse size of ~0.015 pc. The combined radial thickness 
of these sheets or sub-structures in each absorbing region is 
just Nh /riH ~ 10 11 cm if they are at 7? ~ 5 pc, or ~2 x 10 15 
cm at R ~ 1100 pc. 

Some early studies of BALs presented similar ar- 
guments for the existence of dense sub-structures with 
an overall small volume filling factor in quasar outflows 
(|Wevmann et al.lll985l ; lTurnsheldll995l ). Those studies also 
recognized that the creation and survival of such small sub- 
structures pr esent a serious theoretical challenge for out- 
flow models |de Koo 3 Il997h . These concerns about sub- 
structure effective l y end ed when iMurrav et al.l l|l995h and 
iMurrav fc Chianel |l997t) proposed an attractive alternative, 
namely, that the flows could be spatially smooth and con- 
tinuous if there is a shielding medium at the base of the flow 
that blocks most of the quasar's intense X-ray and far-UV 
radiation. With this shielding medium in place, the outflow 
gas behind the shield should be able to maintain the ob- 
served moderate degrees of ionization at much lower gas den- 
sities, thus allowing the flow to have an overall smooth spa- 
tial distribution. Subsequent X-ray observations supported 
this scenario by sh owing that BAL quasa r s are heavily 



absorbed in X-rays (iGreen fc Mathurj Il996l ; iMathur et al.l 



2000; I Gallagher et al 



2002, 2006). 



However, J2123-0050, like other NAL and mini-BAL 
outflow quasars, does not have strong X-ray absorption 
(§1 and §2.2) and therefore it does not have a signif- 
icant X-ray/far-UV sh ield (see also the calculations in 
lHamann fc Simon|[2010l 'l. Small dense sub-structures do ap- 
pear necessary to explain the moderate degree of ioniza- 
tion at the inferred flow locations. Thus the debate about 
the existence of small sub-structures in quasar outflows re- 
mains open. These sub-structures might resemble the small- 
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scale clumps identified recently in one deta iled observationa l 
study of partial covering in a B AL outflow |Hall et al.ll2007f ) . 



4.5 Acceleration & Radiative Shielding: 
Comparisons to BAL Outflows 

The most remarkable aspect of the NAL outflow in J2123- 
0050 is the degree to which it contrasts with BAL flows. For 
example, it has high speeds and ionizations very similar to 
BALs, but the NALs we measure are about a hundred times 
narrower and a hundred times weaker (in REW) than typical 
BAL s (c.f., Table 2 and lWevmann et al.lll99ll ; lKorista et al.l 
I1993T ). Our estimate of the total column density, Nh ~ 
5 x 10 19 cm -2 , is also 1-2 orders of magnitu de less than the 
best available est i mates obtained from BALs (|Hamannl 19981 ; 
lArav et al.ll200ll ; lHamann et~al]|2002t). The weak or absent 
X-ray absorption in J2123-0050 (|just et al.ll2007l) is typical 
of ot h er NAL and mini-BAL outflow sources ( Gibson et al.l 
l20ld ; IChartas et ail |2009| ; iMisawa et ail 120081 1 . but it too 
is very different from the strong X-ray absorption found 
in BA L quasars (jMathur et al 1 120001 : iGallagher et~aLi r2002. 
I2006T ). 

This disparity in the X-ray properties between BAL, 
mini-BAL and NAL outflows has important implications for 
outflow physics. As noted in §4.4, the prevailing models of 
radiatively-driven BAL outflows have a shielding medium 
at the base of the flow that can absorb the quasar's ioniz- 
ing far-UV/X-ray flux and thus preven t the over-ionization 
of the outflow gas beh ind the shield (|Murrav et al.lll995l ; 
iMurrav k, Chiana [l9 97). Moderate ionizations in the out- 
flow gas are essential to maintain significant opacities and 
facilitate the radiative acceleration to high speeds. How- 
ever, the absence of strong X-ray absorption in J2123-0050 
shows that high speeds and moderate ionizations can occur 
in quasar outflows without significant radiative shielding. 

One way to reconcile these results is by a unified scheme 
where the outflows measured via BALs, NALs and mini- 
BALs coe xist at different spatial locations. In situ ations like 
this fe.g.. lGangulv et al.ll200ll ; IChartas et al.ll2009l ). the ac- 
celeration of the NAL and mini-BAL gas might occur in the 
shielded BAL environment (near the accretion disk plane) 
before it moves to its observed location (farther above the 
disk) where the shielding is negligible. The discrete appear- 
ance of NALs and mini-BALs in quasra spectra might be due 
to their formation in discrete blobs or filaments that are cre- 
ated via instabilities along the ragged edge of the ma in BAL 
outflow (|Gangulv et all l200ll ; lHamann et"al] l200Sft . How- 
ever, it is not clear there can be enough vertical force to 
push these blobs significantly away from the disk plane and 
out of their natural radial trajectories away from the UV 
continuum source. A more basic problem with this picture 
is that NALs and mini-BALs have degrees of ionization sim- 
ilar to BALs even though they are not behind the radiative 
shield that is alleged to be important for controlling the 
outflow ionization. If the NAL and mini-BAL gas is acceler- 
ated in a shielded BAL environment where it has moderate 
BAL-like ionizations, how does it maintain those moderate 
BAL-like ionizations after it emerges from the shielded BAL 
zone? 

The obvious answer from our analysis of J2123-0050 
(§3.6 and §4.4) is that the outflow ionization is regulated 
mostly by locally high gas densities and not by a radiative 



shield. High densities could keep the ionization low enough 
for radiative acceleration even if there is no shielding at all. 
Perhaps this is the case generally in quasar outflows, in- 
cluding the BALs. However, high densities also lead to the 
uncomfortable conclusion that the outflows are composed of 
small sub-structures that have an overall small volume fill- 
ing factor (§4.4). The physical processes that might create or 
main tain these sub-structures are not understood |de Kool 
Il997l ). Nonetheless, the evidence for their existence in J2123- 
0050 is at least consistent with high densities playing a more 
important role in the acceleration physics than radiative 
shielding. 



4.6 Energetics & Feedback 

Here we estimate the total mass and kinetic energy in the 
J2123-0050 outflow to determine if it can plausibly be im- 
portant for feedback to the host galaxy's evolution (§1). If 
the flow geometry is approximately like part of a thin spher- 
ical shell, then its total mass is given by 



M 



5 x 10 19 cm- 



R 
5pc 



M . 



(4) 



where Nh ~ 5 x 10 19 cm -2 is our best guess at the total 
column density in all five absorbing systems (§3.6), R ~ 5 
pc is a likely radial distance (§4.1 and §4.4), and Q is the 
global covering fraction of the outflow, i.e., the fraction of 
4-7T steradians cover ed by the flow as seen from the central 
continuum source (Hamann 2000). The value of Q is not 
constrained by our data, but Q ~ 10% is a reasonable guess 
based on estimates of the detection rates of high-velocity 
NAL outflows (§4.2). The outflow kinetic energy defined by 
K = Mv 2 /2 is therefore 



K 



2 x 10 5 



M 



14 Mr. 



12,000 km/s 



ergs 



(5) 



where v ~ 12, 000 km s~ is roughly representative of all 
five systems. 

To estimate the mass loss rate and kinetic energy lu- 
minosity, we would need to make additional ass umptions 
abou t the duration and radial structure of the flow ()Hamannl 
2000). However, if we simply divide K from Equation 5 
by a characteristic flow time, tfi ow ~ R/v ~ 410 yr (for 
R ~ 5 pc and v ~ 12,000 km s _1 ), we obtain a time- 
averaged kinetic luminosity of ~ 2 x 10 42 ergs s _1 
and a ratio relative the quasar's bolometric photon lumi- 
nosity of only (L K ) /L ~ 2 x 10~ 6 (§2.1). This estimate is 
much less than the ratio of Lk/L ~ 5% expected if the 
outflow is to be import ant for feedback (IScannapieco & Oh 
2004 IdT Matte o et al ] 120051 : IProchaska fc Hennawil l2009l l. 
We also note that the total kinetic energy given by Equa- 
tion 5 is equivalent to the non-neutrino output of only ~20 
type II supernovae (at K ~ 10 51 ergs each). 

The total mass and kinetic energy could be much larger 
if the outflow gas is near the maximum radius, R ~ 1100 
pc, derived in §3.6. In that case, Equations 4 and 5 indicate 
M~7x 10 5 M Q and K ~ 10 57 ergs, respectively. However, 
this value of K is still far less than the binding energy of gas 
in a massive galaxy (~10 59 -10 60 ergs), and the increased 
flow time at R ~ 1100 pc leads to a time-averaged kinetic 
energy yield, (Lk) /L ~ 4 x 10~ 4 , that is still two orders of 
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magnitude below the fiducial mark of Lk/L ~ 5% needed 
for feedback. Therefore, at any of the plausible absorber 
locations, the kinetic energy yield from this NAL outflow 
is much too small to be important for feedback to the host 
galaxy's evolution. 



5 SUMMARY 

We discuss multi-epoch spectroscopic observations of five 
outflow NAL systems in the luminous quasar J2123-0050. 
The lines have velocity shifts from —9710 to —14,050 km 
s" 1 and Civ line widths in the range FWHM ~ 62 to 164 
km s _1 (as measured in 2006.64, Table 2). These are the 
highest velocity NALs reported so far to have their forma- 
tion in a quasar outflow confirmed by all three indicators 
of line variability, partial covering of the quasar continuum 
source, and smooth super-thermal line profiles. The five dis- 
tinct absorption line systems require five outflow structures 
with similar physical conditions, similar sizes (covering frac- 
tions) and roughly similar kinematics. All five systems have 
stronger absorption in O VI than C IV, while lower ions such 
as Si IV, Cm and Cn are significantly absent. The observed 
NAL variabilities were well coordinated between the five 
systems, with changes in the line strengths accompanied 
by nearly commensurate changes in the line-of-sight cover- 
ing fractions. The time scale for significant NAL changes is 
^0.63 yr in the quasar rest frame. Our analysis of these data 
provides the following additional constraints on the outflow 
properties: 

1) Absorbing Region Sizes. Partial covering of the UV 
continuum source indicates that the C IV and O VI absorb- 
ing regions have characteristic sizes of order ~0. 01-0. 02 pc 
(§3.4). 

2) Ionization. The absence of C ill lines and stronger 
absorption in O VI compared to C IV in all five systems in- 
dicates an ionization parameter of log U > —0.8 (§3.6). The 
most likely value is near log U ~ —0.4 based on the different 
variabilities observed in Civ and Ovi (§4.3). 

3) Variability. The coordinated nature of the NAL vari- 
ations occurring in ^0.63 yr is best explained by global 
ionization changes in the outflow caused by changes in the 
quasar's ionizing flux (§4.3). The magnitude of the ioniza- 
tion changes corresponds to roughly 0.2-0.3 dex in U. 

4) Metallicity. The secure detection of Lya in one sys- 
tem indicates a metal abundance of [C/H] ~ 0.3±0.2, with a 
lower limit independent of ionization uncertainties of [C/H] 
> 0.2 (§3.7). These results are consistent with other esti- 
mates of quasar metallicities and they support the claim 
that the variable NALs in J2123-0050 form in a quasar out- 
flow. 

5) Total Column Density. We estimate the total hydro- 
gen column density in system C to be Nh ~ 10 19 cm -2 
based on log U ~ —0.4 and the strength of Lya in that 
system. From this and the many similarities between the 
variable systems, we infer that the total column density in 
all five systems is very roughly Nh ~ 5 x 10 19 cm -2 (§3.6). 

6) Dynamics. An upper limit on the radial accelera- 
tion, <3 km s _1 yr - 1 , and an apparent double line-lock in 
the C IV systems suggest that i) the outflow was radiatively 
accelerated, ii) it is now gravitationally unbound and coast- 
ing freely beyond the main acceleration zone, and iii) its 



trajectory is within ~16° of the purely radial (line-of-sight) 
direction (§3.5). 

7) Location. The acceleration upper limit suggests that 
the outflow lines form in a low gravity environment >5 pc 
from the central quasar (§3.5). The minimum gas density 
needed for recombination in ^0.65 yr sets an upper limit on 
the distance of ~1100 pc (in photoionization equilibrium, 
§3.6). Thus the full range of plausible distances is 5 < R < 
1100 pc. If the absorbing structures are blobs or filaments 
created in the inner flow, their short survival times against 
dissipation (in the absence of a confining pressure) indicates 
that they are near the minimum radius of ~5 pc (§4.4). 

8) Density & Sub-Structure. The gas densities corre- 
sponding to distances 5 < R < 1100 pc are 10 8 > ng > 5000 
cm' 3 for a gas in photoionization equilibrium at log U > 
—0.8 (§3.6). These densities are much larger than the aver- 
age value (riff) ~ 215 cm -3 expected from the total column 
density and characteristic size. Therefore, each of the five ab- 
sorbers is composed of dense thin sheets or sub-structures 
with an overall small volume filling factor (§4.4). 

9) Radiative Shielding & Acceleratio n. The absence o f 
strong X-ray absorption in J2123-0050 (|just et al.l 120071 ) 
is typical of NAL and mini-BAL outflows but contrasts 
markedly with BAL quasars. It implies that radiative shield- 
ing in the X-ray/far- UV is not needed to keep the outflow 
ionizations moderate (and BAL-like) and therefore, perhaps, 
it is not needed to facilitate the radiative acceleration of the 
flow to high speeds. We argue that the ionization is con- 
trolled, instead, by locally high gas densities in small outflow 
sub-structures (§4.5). 

10) Energetics & Feedback. At the most likely radial 
distance near R ~ 5 pc, this outflow has a total mass of 
only ~14 M© and kinetic energy ~2 x 10 52 ergs. At any 
of the plausible larger distances noted above, the kinetic 
energy yield from this outflow is still at least two orders of 
magnitude too small to be important for feedback to the 
host galaxy's evolution (§4.6). 
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APPENDIX A: PHOTOIONIZATION 
CALCULATIONS 

We use th e numerical code C loudy, version 08.00 (last de- 
scribed bv lFerland et al.l Il998l ') . to compute ionization frac- 
tions for clouds that are in photoionization equilibrium with 
a quasar continuum source. The results have wide applica- 
bility to quasar outflows and a variety of other near-quasar 
environments. The calculations assume the gas clouds are 
dust free, at a constant density throughout, optically thin in 
the Lyman continuum and geometrically thin such that their 
radial thickness, AR, is small compared to their distance, R, 
from the ionizing emission source, i.e., AR/R <C 1. The op- 
tically thin assumption ensures that the model clouds have a 
uniform ionization structure throughout. This is appropriate 
for the majority of quasar outflow systems, including J2123- 
0050, that are highly ionized with AT (Hi) < 10 172 cm -2 . 
(Moderate continuum opacities can occur at the He II ioniza- 
tion edge for 10 16 ' 5 < AT (Hi) < 10 17 ' 2 cm" 2 , but this has no 
significant affe ct on the ionization fractions discussed here, 
lHamanr]|l997h The calculations assume that the metallicity 
is solar and the total hydrogen density is uh = 10 5 cm -3 , 
although the specific v alues of these p arameters do not affect 
the ionization results (|HamanrJll997f ). 

The quasar spectrum is a simple segmented powerlaw 
consistent with observati onal constraints across the most im- 
portant far-UV ener g ies jElvis et al.ll 19941; iTelfer et al.ll2002l; 
Steffen et al.l 120061 ; iRichards et all 120061 ; iHopkins et all 



20071 . and refs. therein). It has the form /„ oc v a with slope 
a = -1.0 from 0.5 eV to 13.6 eV, a = -2.56 from 13.6 eV to 
136 eV, and a — — 1.0 from 136 eV to 136 keV. Sharper cut- 
offs are applied to the high and low energy extremes (a = 3 
below 0.5 eV and a = —3 above 136 keV) that have no 
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affect on the computed ionizations. This spectrum yields 
ctox = —1.6 (appropriate for luminous r adio-quiet quasars, 
IStrateva et~al1 boOa : ISteffen et aDl2006h . and a bolometric 
correction factor at 1450 A of L w 4.4 i/Z/„(1450 A) (see also 
IWarner et afll2004D . 

We describe the ionization state of the gas in 
terms of the dimensionless ionization parameter, U = 
&h /(Attcuh R 2 ), where &h is the h ydrogen-ionizing pho - 
ton luminosity of the central source (|Ferland et al l 1998). 
For the continuum shape described above, this yields the 
following expression for the radial distance 

R = 225 1/2 ( iJf^L) 1/2 (M) 1/2 pc (A1) 

v 10 47 ergs a- 1 J \ n H J \Uj y ' 

Figure Al shows the calculated ionization fractions for 
Hi and several metal ions commonly observed in near-UV 
spectra of quasars, namely, ions of C, N, O and Si. One im- 
portant application of these results to absorption line work 
is derivations of metal-to- hydrogen (M/H) abundance ratios 
based on this relation 

where (H/M)q is the solar abundance ratio by number, and 
N and / are the column densities and ionization fractions, 
respectively, for Hi and the metal M in ion stage i. Ideally, 
one would apply Equation 2 after comparing the measured 
column densities in several metal ions to the calculated re- 
sults (Figure Al) in order to define U and an appropriate 
ionization correction factor, /(Hi) //(Mi). However, if ion- 
ization constraints from the data are not available or not 
reliable, one can also use minimum values of the ioniza- 
tion corrections in Equation A2 to derive fi rm lower limits 
on the metal-to-hydrogen abundance ratios |Hamannll 19971 ; 
lHamann fc Ferland|[l999l 'l. 

Table Al lists the minimum correction factors for vari- 
ous metal ions, log(/(Hi)//(Mj)), together with the values 
of the ionization parameter, log U, where these minima oc- 
cur. This information is provided in the table according to 
the element (column 1) and ionization stage II- VI (columns 
2-6) for ions that are sometimes detected in quasar outflow 
spectra. Note that the minimum correction factors occur 
near, but at larger U than, the peaks in the /(Mi) curves 
shown in Figure Al. For convenience, we note that the most 
recent solar abundance ratios are log(H/M)e = 3.57 for 
carbon, 4.17 for nitrogen, 3.31 for oxygen, 5.55 for alu- 
minum , 4.49 for silicon, 6. 59 for phosphorus and 4.88 for 
sulfur l|Asplund et alj|2009h . 
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Figure Al. Ion fractions, /(Mi), for selected stages of the ele- 
ments C, N, O and Si plotted against the ionization parameter, 
logf/, in optically thin photoionized clouds (see text). From left 
to right the three curves plotted for carbon are for C II, C III and 
Civ. The remaining curves (left to right) are for nitrogen: Nm, 
Niv and Nv, oxygen: Om, Oiv, Ov and Ovi, and silicon: Sin, 
Sim and Siiv. The Hi fraction, log /(Hi), at each logU is marked 
along the horizontal bar near the top. 



Table Al. Minimum Ionization Corrections 



Element 


II 


III 


IV 


V 


VI 




minimum log(/(H 


0//(Mi)): 






Carbon 


-1.75 


-3.32 


-3.80 


-5.39 






Nitrogen 


-1.75 


-3.28 


-4.15 


-4.55 


-5 


89 


Oxygen 


-1.81 


-3.02 


-4.17 


-4.69 


-5 


06 


Aluminum 


-2.12 


-2.20 


-3.36 


-4.11 


-4 


96 


Silicon 


-1.65 


-2.93 


-3.03 


-3.89 


-4 


51 


Phosphorus 


-1.90 


-3.10 


-3.34 


-3.57 


-4 


12 


Sulfur 


-1.50 


-2.95 


-3.73 


-3.73 


-3 


97 




log U at 


minimum 


log(/(H 


l)//(Mi)) 






Carbon 


-2.60 


-1.55 


-0.75 


0.30 






Nitrogen 


-3.20 


-1.50 


-0.75 


-0.15 





65 


Oxygen 


-3.15 


-1.65 


-0.75 


-0.15 





30 


Aluminum 


-2.60 


-2.10 


-1.50 


-0.70 


-0 


10 


Silicon 


-2.80 


-1.95 


-1.55 


-1.00 


-0 


35 


Phosphorus 


-2.60 


-1.80 


-1.35 


-1.00 


-0 


50 


Sulfur 


-2.90 


-1.65 


-1.20 


-0.85 


-0 


60 
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